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ABSTRACT  This  article  studies  a  type  of  new  model  strap 
down  terminal  guidance  system.  It  contains  a  strap  down  type 
guidance  head,  a  set  of  strap  down  type  inertial  sensors,  and  in 
the  three  speed  gyroscope  and  three  line  accelerometer  are 
partially  shared  measurement  components  associated  with 
navigation,  guidance,  and  piloting  instruments.  The  article 
presents  a  strap  down  terminal  guidance  design  associated  with 
the  estimation  of  target  status  relative  to  missile  movements  in 
inertial  coordinate  systems,  deriving  the  structures  and 
mathematical  models  associated  with  the  types  of  systems. in 
question.  It  puts  forward  design  algorithms  for  generalized 
Karman  (phonetic)  filters  as  well  as  optimal  guidance  patterns. 
Finally,  use  is  made  of  digital  simulations  to  empirically  verify 
the  feasibility  of  the  systems  in  question  to  handle  close  range, 
highly  maneuverable  targets . 

KEY  WORDS  Strap  down  terminal  guidance  system,  Strap 

down  type  guidance  head,  Generalized  Karman  (phonetic)  filter. 
Optimized  guidance. 


SYSTEM  STRUCTURE  AND  OPERATIONAL  PRINCIPLES 

In  order  to  be  appropriate  to  the  requirements  associated 
with  countering  close  range  (under  5  kilometers)  highly 
maneuverable  (normal  overload  reaching  9g)  aerial  targets,  on  the 
one  hand,  use  is  made  of  guidance  head  technology  progress 
achieved  in  the  area  of  expanded  instantaneous  fields  of  vision. 
It  is  possible  to  make  guidance  heads  get  rid  of  universal  frames 
and  fix  them  directly  to  missile  bodies.  This  kind  of  strap  down 
type  guidance  head  not  only  makes  tracking  angular  velocities 
almost  unlimited.  It  also  increases  system  reliability.  On  the 
other  hand,  guidance  system  design  bends  every  effort  to  present 
new  structural  arrangements,  collecting,  and,  in . conjunction  with 
that,  organizing  and  extracting  as  much  information  as  possible, 
using  microcomputers  to  introduce  control  systems,  and 
constructing  new  models  of  guidance  patterns  in  order  to  realize 
optimal  guidance.  The  strap  down  terminal  guidance  system  which 
this  article  studies  opts  for  the  use  of  a  strap  down  type 
guidance  head  as  well  as  a  set  of  strap  down  type  inertial 
sensors  (including  three  velocity  gyroscopes  and  three  linear 
accelerometers)  to  act  as  measurement  components.  It  is  a  new 
type  of  design  facing  precisely  toward  future  combat. 

The  system  structural  diagram  can  be  seen  in  Fig.l.  The 
three  gyroscopic  velocity  measurements  associated  with  angular 
velocities  p,  q,  and  r  along  the  x,  y,  and  z  axes  of  missile 
body  coordinate  systems  go  through  directional  cosine  matrix 
integration.  It  is  possible  to  obtain  coordinate  transform 
matrices  between  missile  body  coordinate  systems  and  inertial 
coordinate  systems .  Accelerometers  measure  the  three  linear 
missile  accelerations  AmXb,  Amyb,  and  Amzb  along  the  three  a^es 
of  the  missile  coordinate  system.  Strap  down  type  guidance  hea  s 
measure  the  angles  of  sight  Apb  and  the  azimuth  angles  Ayb  of 
targets  in  missile  coordinate  systems.  In  all  cases,  they  go 
through  coordinate  transforms  to  create  the  linear  accelerations 
Amxi  Amyi  and  Amzi  as  well  as  the  corresponding  angles  Api  and 
Avi  in  inertial  coordinate  systems.  Use  is  made  of  Amxi,  Amyi, 
and  Amzi  to  act  as  input  signals.  Api  and  Ayi  act  as  measuremen 
signals.  Going  through  generalized  Karman  (phonetic)  filters,  it 
is  possible  to  do  estimates— in  inertial  coordinate  systems  - 
of  various  parameter  vectors  such  as  missile  target  relative 
positions  SR  ,  relative  velocity  VR  ,  as  well  as  target 
acceleration  AT  ,  and  so  on.  These  parameters  are  needed  m  the 
formation  of  optimum  guidance  commands.  It  is  still  necessary  to 
take  these  parameters  and  transform  them  into  missile  coordinate 
svstems.  Use  is  made  of  the  x  axis  components  as  well  as 
measured  missile  acceleration  values  along  missile  axes  m  order 
to  construct  a  type  of  relatively  accurate  flight  time  remaining 
tgo  algorithm.  The  tgo  in  question  is  an  important  parameter  in 
the  formation  of  optimum  guidance  commands.  Using  the  missile  _ 
guidance  commands  obtained  in  inertial  coordinate  systems,  Amcxi, 
Amcyi,  and  Amczi— -through  coordinate  transformations  formation 
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is  done  of  guidance 


commands  Amcyb  and  Amczb  (missile  axis 


Fig.l  Guidance  System  Composition  Schematic 

*ig  =-1-0  t  i  lleaible)  Transform  (3) 

Transform^*  <^OeKed  ™  phonetic)  ^^tation 
(9)  Coordinate  (illegal  >  Velocity  Gyroscope 

<“>  Direct ional^ Cos ine < Matrix  Integrate 


directions  not  controlled)  in  missile  coordinate  systems.  Then, 
amplitude  limitation  is  applied  to  facilitate  the  output  of  the 
pitch  channel  and  yaw  channel  commands  Ap  and  Ay  to  the  automatic 
pilot  in  order  to  control  missile  flight  toward  the  target. 

MATHEMATICAL  SYSTEM  MODELS 


I.  Coordinate  Systems  Used  and  Their  Transform 
Relationships 

As  far  as  missile  fuselage  coordinate  systems  are . concerned, 
they  are  firmly  connected  to  the  missile  body. _  The  origin  point 
lies  on  the  center  of  mass  of  the  missile,  x  is  the  longitudinal 
axis  along  the  missile.  The  2  axis  lies  perpendicular  to  x  m  a 
plane  of  symmetry  longitudinal  along  the  missile.  See  Fig. 2. 


Fig. 2  Missile  Body  Coordinate  System  as  well  as  Angles  of  Sight 
and  Azimuth  Angles 


As  far  as  inertial  coordinate  systems  are  concerned,  they 
are  ftLd  Movable  relative  to  the  ground.  They  are  a  reference 
system  associated  with  the  movements  of  missiles  and  targets. 
lire option  is  made  for  the  missile  body  coordinate  system  at 
the  instant  of  launch  to  act  as  inertial  coordinate  system. 

Setting  the  coordinate  transform  matrix  T  s  from  mer  1a 
coordinate  system  to  missile  body  coordinate  system  at  the 
instant  t  to  be  (  £ t  / 2  £  , 
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the  three  lines  of  1  stand  respectively  for  the  directional 
cosines  associated  with  the  three  inertial  coordinate  system 
axes  in  the  missile  body  coordinate  system.  Therefore,  they 
have  the  designation  of  direction  cosine  matrix. .  Giving 
consideration  to  microrotational  movements  associated  with  the 
missile  body,  it  is  possible  to  deduce  the  microequations  which 
should  be  satisfied  by  direction  cosine  matrices 
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Tn  the  eauations,  p,  q,  and  r  are  measured  by  velocity  . 

gyroscopes .  MaKing'ule  of  fourth  order: (phonetro) 
methods  to  make  solutions,  it  is  possible  to  obtain  (T , 
any  instant. 

II.  Strap  Down  Type  Guidance  Heads 

As  far  as  measuring  angles  of'  sight  and  position  angles  of 
taraets  in  missile  body  coordinate  systems  is  concerned,  see 
TPin^?  inb  is  the  included  angle  between  the  projection  of 
.  “[*  (3t-_Tnissil e  relative  position  vector  RTM  in  the  plane  xoz  and 

the^missile^bod^longitudinal  axis  x.  Xyb  is  the  included  angle 
between  ihe  projection  of  RTM  in  the  plane  xoy  and  x.  By 
geometrical  relationships,  one  has 
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*  Ryb 

xyfc-ti  1~T+Vyi> 


In  equations,  Rxb,  Ryb, 
of  RTM  in  missile  body 
measured  noise. 


and  Rzb  are  the  various  axial  projections 
coordinate  systems.  Vpb  and  Vyb  are 


III.  Conversion  of  Guidance  Head  and  Linear  Accelerometer 
Information 

rnmnlptelv  or  making  use  of  Karman  (phonetic) filters  to 
calculate^estimates  for  mfssile-target  relative  motion  states  in 
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inertial  coordinate  systems,  it  is  necessary  to  take  information 
measured  by  guidance  heads  and  linear  accelerometers  in  inertial 
coordinate  systems  and  convert  it  into  information  in  inertial 
coordinate  systems .  Due  to  the  orthogonal  nature  of  direction  _ 
cosine  matrices,  the  coordinate  transfer  matrix  X,  from  missile 
body  coordinate  systems  to  inertial  coordinate  systems  is  nothing 

else  than 


Therefore, 
systems— -in 
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the  vectors  Amb  and  RTMb  in  missile  body  coordinate 
inertial  coordinate  systems^-are 


Am  j.  ~  1  i  *  AM<b 

RjMi 


(5) 


The  three  component  vectors  associated  with  RTMi  and  RTMb  are 

Rnds^  Rxi  %i 

RrMbsC3n  %b 

Therefore,  one  has 

Rxi=^*Rxb  +  in>Ryb+n!Rz'b  ') 

Ry  j_  =  £  i  ^xb  +  m  2  %  b  +  n  2'  Rz  b 
H2  ~£ s^xb  +  m  5Ryb  *n  sRz  b  • 
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Because  of  this,  angles  of 
guidance  heads  in  inertial 


sight  and  position  angles  measured  by 
coordinate  systems  are 
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By  the  same  reasoning, 
respectively 


the  three  components 


of  AMi  and  AMb  are, 


AMi-tA«xi  A“Y  i  A  M  z  i  ^  * 
AMbs£A“xfc  A«l*fc 


(9) 


By  contrast,  the  v^ous  co^^associ^^  Unear 

accelerometers  m  inertia 
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IV.  Besigns  Associated  with  Generalised  Kar^an  (phonetic, 
Filters  ^  the  of  optimum  gnidance,  there 
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is  a  need  for  information  associated  with  all  the  various  vectors 
for  missile-target  relative  positions  SR  ( SR  =  RTM) ,  relative 
velocities  VR  ,  as  well  as  target  accelerations  AT  .  Moreover, 
these  status  variables  are  not  able  to  be  measured  by  systems  on 
missiles.  There  is  a  need  for  additional  calculations.  Use  is 
made  of  strap  down  guidance  head  linearization  as  well  as 
information  measured  by  linear  accelerometers  in. inertial 
coordinate  systems .  Moreover,  giving  consideration  to  the 
existence  of  noise,  there  is  a  need  to  design  generalized  Karman 
(phonetic)  filters  in  order  to  calculate  these  states  in  inertial 

coordinate  systems.  . 

The  micro  equations  which  the  three  vectors  SR  ,  VR  ,  and  AT 
which  await  solution,  satisfy  are  determined  by  missile-target 
relative  motion  equations.  That  is,  equations  of  state 
associated  with  generalized  Kanru'n  (phonetic)  filter  problems  are 


Ss  =  V, 

V *  =  -Au  *  I,,, 

a 

Ay  =  ~  A.  A  r  t  i,? 


(ID 


In  the  equations,  WM  and  WT  are,  respectively,  acceleration 
SoceSs  noises  associated  with  missiles  and  targets.  Here 
described  target  mechanical  characteristics  are  shown  using  firs 
order  Maerkefu  (phonetic)  processes,  written  as  matrices,  one 

forms 


X  =  r  T  0  A  k; 


(12) 


In  equations, 
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In  order  to 


facilitate  digital  control,  the  particularized  state 
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equation  is 


xck  +  i  )  -0  ck+ 1„  kJxdO+Gck+l,  k  )Am  ck  > 
+rck+l,  lOWcki 


(13) 
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In  equations, 


'US  ,  — ^  o  ^  t 
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From  the  point  of  view  P*  ™*™^*'  Record  it 
x  and  expand  it.  In  conjunction  with  that, 


x  s  £  S  a  Y  a  -A-r  3  r  s  £  S  s  x^  s  iyi  R  r  z  i  ViXi  "^syi 
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=  £x  j  x s  x  3  —  x , 3  r 


Moreover , 


£SBXi  Ssyi  Sszi3  =CHxi  Ryi  Kzi3 


Because  of  this,  strap  down  type  guidance  head  equations  can  be 
written  into  the  matrix  form 

(16) 

2^hcx3-s-Y 
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In  equations, 
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For  the  sake  of  its  linearization  and  particularization,  around 


A 

predicted  state  + 

series,  and,  in  conjunction  with  that, 
and  higher  order  terms.  It  is  possible 
observational  equation 


expand  a  Taylor 
omit  the  various  second 
to  obtain  the 


2Ck  +  I>  -Hck+Dxck  +  D+Y  ck  +  l)+vck+  I  ) 


In  equations,  the  observational  matrix  H(k+1)  is 
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After  obtaining  linearized  and  particularized  equations  of 
state  and  observational  equations,  on  the  basis  of  the  theory  of 
optimization  calculations,  it  is  possible  to  obtain  the 
generalized  Karman  (phonetic)  filter  equations  possessing  input 
signals  (AM  (k) )  as  the  equation  set  presented  below 


xcK*2S k)  «0ck-a-2o  lOxck/k) 
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(18) 


In  equations 

Qjj  s  EC  Wjj  *  w£  3  s  CVj£  «  V-^3 
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V.  Derivation  of  Optimum  Guidance  Patterns 


In  order  to  facilitate  initial  research — when  deriving 
optimum  guidance  patterns — automatic  pilots  and  missile  body 
inertias  are  temporarily  ignored.  Because  of  this,  equations  of 
state  associated  with  controlled  objects  are  then  nothing  else 
than  target-missile  relative  motion  equations.  Taking  the  vector 
microequations  in  question  and  expanding  them,  it  is  possible  to 
obtain  9  scalar  equations.  They  are  capable — in  accordance  with 
the  three  axes  of  inertial  coordinate  systems — of  being  divided 
into  three  microequation  sets  associated  with  solution  pairs.  In 
all  cases,  each  set  of  equations  has  the  general  form  /5-7 
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Selection  is  made  of  the  second  order  form  of  performance  index 
function  set  out  below 
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In  equations, 
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It  is  possible  to  solve  for  optimum  control  as 
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That  is  nothing 
along  the  three 


else  than  to  say  that  optimum  gtidanoe  commands 
inertial  coordinate  system  axes  are 
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in  order  to  make  target  n.iss^uantitles^as^mall^as.possrble^ 

|^etChi°useSh:£1opt^  guidance  patterns  which  are 
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similar  guidance  patterns 
which  are 
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missile  body  coordinate 


systems , 
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• _  y  it  is  fi  one 

Setting  out  ^.“it  is  possible  to  solve 

SeSe°“icSation  formula  for  tgo  as 
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In  equations, 
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transformations  into  corresponding  components  in  missile  axis 
directions .  Again  making  use  of  acceleration  values  measured  in 
missile  axis  directions,  it  is  then  possible  to  calculate  tgo  . 
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VII.  Formation  of  Control  Commands 

Control  commands  obtained  in  inertial  coordinate  systems 

need  to  go  through  coordinate  transfer  matrix  j*a  to  be  control 

commands  in  missile  body  coordinate  systems.  In  conjunction  with 
this,  consider  the  missile’s  largest  normal  acceleration  to  be  A. 
Then,  amplitude  limited  control  commands  outputted  to  automatic 
pilots  are 


Ap  -  A 


\ 


Ay  ®  A 


C  I 


(26) 
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NUMERICAL  SIMULATION  RESULTS 


In  order  to  study  the  performance  of  the  strap  down . terminal 
guidance  design  in  question,  ballistic  numerical  simulation 
calculations  were  carried  out  on  a  few  typical  target  motion 
situations.  In  simulations,  it  is  assumed  that  automatic  pilot 
and  missile  body  inertias  are  not  calculated.  Velocity  gyroscope 
outputs  are  regarded  as  a  white  noise  process.  During 
simulations,  following  along  with  alterations  in  missile-target 
relative  distances  r,  guidance  system  sampling  cycles  At  are 
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On  the  basis  of  reference  [1],  the  power  spectrum  density 
matrices  R(t)  and  Q(t)  associated  with  observation  noise  and 
process  noise  in  generalized  Karmen  (phonetic)  filters  are 
selected  to  be  ' 
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Initial  values  for  relative  missile-target  positions  as  well 
as  relative  velocities  are  replaced  with  actual  values . _  Target 
velocity  initial  values  are  0  (illegible),  X  -  1.  The  initial 
value  P(%)  of  calculation  error  variance  matrices  is 

P(%)  =  diag  {000111  2500  2500  2500} 

Two  typical  target  movement  situations  are  presented  here. 
Situation  I.  Target  is  in  linear  flight  at  uniform  speed  m 
inertial  coordinate  system.  Velocities  along  x,  y,  and  z  axes 
are,  respectively,  300,  50,  and  50  m/s.  The  launch  distance  is 
2500m.  Situation  II.  Target  is  in  high  maneuver  flight  within 
inertial  coordinate  system.  Accelerations  along  y  and  z  axes 


are  both  50m/s2  .  Initial  velocity  is  300  m/s.  In  the  two  types 
of  situations,  missiles  are  aimed  and  launched  along  the  x  axis. 

MlSS1Targetemiss  quantity  calculation  results  in  inertial 
coordinate  system  are 
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(1)  Situation  (2)  Overall  Amount  (3)  Components 
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Fig. 3  presents  estimated  value  curves  associated  with 
Situation  I  generalized  Karmen  (phonetic)  filters.  In  order  to 
facilitate  comparison,  actual  values  are  drawn  on  the  same  graph. 
Situations  of  change  associated  with  missile  control  commands  are 
also  given.  Fig. 4  is  calculation  results  corresponding  to 

Situation  II.  . 

Conclusions.  The  strap  down  terminal  guidance  design  m 
question  is  appropriate  for  countering  close  range,  highly 
maneuverable  targets.  The  digital  models  are  reasonable. 
Generalized  Karmen  (phonetic)  filters  are  capable  of  relatively 
accurately  estimating  missile-target  relative  movement  states  as 
well  as  target  movement  states.  Methods  for  estimating  time 
remaining  are  comparatively  good.  Optimum  guidance  patterns  are 
workable.  With  regard  to  optimum  guidance  problems  associated 
with  the  consideration  of  automatic  pilot  and  missile  bo  y 
inertias,  it  is  planned  to  carry  out  research  m  another  article. 


Fig. 3  (a)  Key:  (1)  Actual  Value  (2)  Estimated  Value 


Fig. 3  (b)  Key:  (1)  Actual  Value  (2)  Estimated  Value 
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Fig. 4(b)  (1) 


Actual  Value  (2)  Estimated  Value 


Fig. 4(c)  (1) 


Actual  Value  (2)  Estimated  Value 
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